Wigner molecules in quantum dots 
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We perform unrestricted Hartree-Fock (HF) calculations 
for electrons in a parabolic quantum dot at zero magnetic 
field. The crossover from Fermi liquid to Wigner molecule 
behavior is studied for up to eight electrons and various spin 
components S^. We compare the results with numerically 
exact path-integral Monte Carlo simulations and earlier HF 
studies. Even in the strongly correlated regime the symme- 
try breaking HF solutions provide accurate estimates for the 
energies and describe the one-particle densities qualitatively. 
However, the HF approximation favors the formation of a 
Wigner molecule and produces azimuthal modulations of the 
density for even numbers of electrons in one spatial shell. 

PACS numbers: 73.20. Dx, 71.10Ay, 71.10.Hf 



The last decade has seen an enormous interest in quan- 
tum dots, i.e., a small number of 2D electrons confined 
in a semiconductor heterostructure Experimentally, 
the A^-electron states of these systems are studied by 
means of far- infrared capacitance |^ and transport 
spectroscopy [Q, and exhibit features of quantization of 
charge and energy. Theoretically, a whole arsenal of 
methods for interacting electronic systems together with 
increasing computational power applies: Exact diagonal- 
ization techniques ||], density functional theory |^ and 
quantum Monte Carlo methods j7|. Recently, also the 
very strongly correlated regime of small electronic densi- 
ties has attracted considerable interest [p|-p^ . In particu- 
lar, the formation of a Wigner molecule was studied with 
Quantum Monte Carlo (QMC)J3, Hartree-Fock (HF) ||] 
and configuration interaction mQ\ calculations. 

In this paper we reconsider the Hartree-Fock approxi- 
mation focusing on the crossover to the Wigner regime. 
An important practical and conceptual question to be 
clarified is how correlations beyond the mean field ap- 
proximation contribute to the exact energy and whether 
the Wigner molecule is described reliably within HF. In 
its unrestricted version, allowing for symmetry-broken 
solutions, HF approaches the true ground state energy 
considerably better than restricted HF which preserves 
the rotational symmetry of the Hamiltonian. This is 
achieved, however, on expense of the quality of the wave 
functions. 

In the case of a strong central potential in three di- 
mensions, such as in real atoms, the HF approximation 
is known to yield useful results both for energies and 
wave functions. Here, we show by comparison with exact 
Monte Carlo data M that even for strong interaction un- 



restricted HF can give very good estimates for the ground 
state energies of 2D quantum dots. Only tiny energy dif- 
ferences between different spin states cannot be resolved 
reliably. On the other hand, the charge density distribu- 
tion resulting from unrestricted HF cannot quantitatively 
describe the strongly correlated regime. For even num- 
ber of electrons per shell, the HF densities show effects of 
localization due to the strong electron interaction. Fur- 
thermore, within HF |^ this crystallization sets in rather 
too early at higher densities than in the QMC study [||. 

Model — We study a two-dimensional parabolic quan- 
tum dot with N electrons at zero magnetic field as it has 
been discussed by many authors [5-10]. Measuring en- 
ergy in units of the oscillator energy hujo and length in 
units of Iq — ^JhJmFujQ, where m* is the effective mass, 
the dimensionless Hamiltonian reads 
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Here we have introduced the dimensionless coupling con- 
stant A ~ Iq/ciq — e^ / kIqUuiq with the effective Bohr 
radius ag and the dielectric constant n. For example 
A = 2 corresponds to Tiuj^ ~ 3meV for a GaAs quantum 
dot. Since this Hamiltonian is rotationally invariant and 
spin independent, the exact eigenfunctions can be cho- 
sen as simultaneous eigenfunctions of the total angular 
momentum the total spin S'^^^ and its z-component 
5*°*. These eigenfunctions and corresponding densities 
are then rotationally invariant. 

Method — HF theory consists in approximating the 
many-particle wave function by an optimal single Slater 
determinant 
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This HF Slater determinant is build up by single par- 
ticle orbitals ^Pi{r), which are expanded in the angular 
momentum basis of the 2D harmonic oscillator (Fock- 
Darwin states) 
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Here, n and M are the radial and angular quantum num- 
bers, and >c[f^' is a Laguerre polynomial. In the unre- 
stricted HF approximation an orbital has the expansion 

ip,{r) = <M{r\nMa,) , (4) 

71=0, OC 

M = — oo , oc 

where Ui = ± ^ is the fixed electron spin of the i-th or- 
bital. This orbital is no longer an eigenfunction of the one 
particle angular momentum, therefore the HF Slater de- 
terminant is in general not an eigenfunction of , but 
a deformed, symmetry broken solution Similarly, 
the HF solution is in general not an eigenfunction of the 
total spin but only of its z-component with eigenvalue 

Minimizing the HF energy E^^ ^ (*"^| H j^f^^) and 
imposing orthonormality between the HF orbitals yields 
the HF equations for the expansion coefficients u\^j^,i 

Y.{{l\h \a) + ^(7/?k|a/3')P/3'/3}< = • (5) 

a 13(3' 

Here, greek indices abbreviate the quantum numbers of 
(^, for example a = {n,M,a). The density matrix is 

defined by paa' — 'Ki'^^a* ^ (tI ^ \^) one par- 

ticle matrix element and {aa'\w\(3l3') = {aa'\w\l3l3') — 
{aa'\ w |/3'/3) denotes the antisymmetrized Coulomb ma- 
trix element, which we calculated analytically by trans- 
formation into relative and center of mass coordinates 
p3[ . The nonlinear selfconsistent eigenvalue problem (j^) 
for the expansion coefficients is solved by iterative di- 
agonalization starting from an initial guess for the den- 
sity matrix. Several starting guesses have to be chosen 
in order to avoid local minima. The true minimum can 
also be identified by its one particle density, 

N 

n«F(r)=^|^.(r)P, (6) 

4=1 

which in the strong coupling limit mirrors the geome- 
try of classical electrostatic point-charges as we discuss 
below. 

In our calculations we used up to 55 Fock-Darwin 
states for each spin direction. We present results for var- 
ious sets of N, A, Sz from = 2 to 8 particles covering 
the whole range of interaction strengths from A = 2 to 
10. 

HF energies and spin — The most important proper- 
ties of the HF ground state energies can already be seen 
for N — 2 (quantum dot Helium). Here one knows the 
exact solution from diagonalization of the Hamiltonian 
for the relative motion We have reproduced these 
results in Table || The true ground state is always a 
singlet and the energy gap to the triplet vanishes slowly 
as A — !■ oo. In contrast, the HF approximation finds a 
triplet ground state for A_> 5, and the energies of singlet 



and triplet converge more rapidly |lj]. The (absolute 
and relative) error is largest in the case oi Sz = and 
A 2. In the polarized case the error of HF increases 
monotonically with A. These features of the HF energies 
can be seen for all N. 

For > 2 we compare with the path integral QMC 
data of Ref. , which have been obtained for the very 
low temperature T = O.lhtoo/kB- The HF ground state 
energies are always above the QMC energies for all par- 
ticle numbers investigated so that the QMC data may 
serve effectively as zero temperature reference points. In 
the case of iV = 5 electrons, QMC gives a ground state 
spin of 1/2. In Fig. |^ we plot the energy difference of 
the different spin states with respect to the true ground 
state, AE — Es — E^J!^'^ , for various spin states as a 
function of the coupling strength A. Again we see that 
AE is largest for A « 2 and smallest spin. It remains 
constant for A >, 4, resulting in a very low relative error 
of only 1 — 2% in the low density regime. The HF energy 
is very accurate in the case of S'z = 5/2 which becomes 
the HF ground state for A >, 4. This unphysically high 
ground state spin is due to the exchange term in the HF 
energy, which lowers the energy only for parallel spins. 
Therefore the spin ordering of the HF energies is just 
reversed. 

For A?^ = 8 we give the HF energies in Table ||. The 
QMC predicts a transition of the ground state spin from 
5 = 1 to S* = 2 for A >, 4. Again HF finds the wrong spin 
ordering in the correlated regime. The relative error is 
largest for the unpolarized states. The convergence of the 
HF energies for different 5*2 for A >, 6 can be understood 
in the classical picture of localized electrons without over- 
lap and therefore no spin sensitivity. 

HF densities — Next we study the HF one particle den- 
sities (^), where the crossover to the Wigner molecule 
is discernible. However, one has to keep in mind that 
the deformation and structure of the HF densities arise 
from the symmetry violating mean field and therefore 
may be artificial (whereas the HF energies are a true 
upper bounds for the ground state energies). Second we 
want to point out that a symmetry breaking Slater deter- 
minant does not necessarily mean that its corresponding 
density displays a molecule-like structure. In Fig. ^ we 
display the densities for A^ = 2 to 5 for strong coupling 
A = 6 and maximal spin Sz = N/ 2 (From A >, 6 the 
HF densities are essentially the same for all spins.). In 
the case of 2 and 4 electrons they show quite distinct az- 
imuthal maxima. This is not the case for the odd electron 
numbers 3 and 5: Surprisingly the HF densities, though 
belonging to a deformed Slater determinant, seem to be 
rotationally symmetric. 

In order to understand this even-odd effect consider 
an exact spin-polarized A^-electron wave function '5 n for 
the Wigner molecule case. An even number of electrons 
in one spatial shell carries a nonzero angular momen- 



2 



turn ±hN/2 In this fashion the HF wavefunctions 

with modulated density {N = 2, 4) can be interpreted as 
standing waves from a superposition of opposite angular 
momenta [ [l5[ |. 

However the HF densities display the right classical 
filling for the spatial shells: From the maxima in the 
densities one can also read off the Brueckner parameter 
rs, which is defined as the the nearest-neighbor distance 
in units of Og. These values agree well with those 
obtained by QMC and from a model of classical point- 
charges. In Fig. Ijwe show a 7-electron Wigner molecule 
for two interaction strengths. For A = 2 the HF density 
is already 6-fold modulated with a central electron. Con- 
fronting this with the exact density in Fig. 1 of Ref. 
we notice that the crystallization (as suggested by the HF 
densities) occurs too early. This explains why the authors 
of Ref. observe Wigner crystallization for higher den- 
sities than in the exact QMC study js) . For higher A the 
molecule with two spatial shells becomes more and more 
distinct. The maxima agree very well with the classical 
formula = A** (2.25 -I- 1/a/3) for the 6-fold geometry. 

Finally we show in Fig. || the crystallization for A*" = 6 
and Sz = 0- This unpolarized case depicts a crossover 
from a 6-fold arrangement [Fig. ^(a)] to a 5-fold geom- 
etry [Fig. ^(c) and (d)]. The 6-elcctron molecule was 
also studied by Reimann et al. [ic| ], who found by means 
of configuration interaction calculations that the true 
ground state was unpolarized with a 6-fold symmetry 
up to at least A « 3.5. Within HF the Sz = state 
acquires a 5-fold symmetry already for A > 2.85. In this 
range the HF density is distorted [A = 4 in Fig. |^(c)] and 
then again apparently round for higher lambda [A = 6 in 
Fig. ^(d)] with a central maximum. Fig. ^(b) shows the 
6-fold isomer for A 4, which is by 0.33 higher in energy 
than the ground state in Fig. §(c). In contrast the spin- 
polarized state exhibits 5-fold symmetry throughout the 
whole parameter range [ll6|| . 

In conclusion we have shown in how far an unrestricted 
HF description of the Wigner molecule in quantum dots 
is reliable. The HF energies are good estimates of the 
true ground state energies especially for spin-polarized 
states. The energy differences for different spin-states in 
the strong interacting regime cannot be resolved prop- 
erly. We find deformed HF densities in the regime of 
intermediate interaction up to A <, 4. For strong cor- 
relation the densities are azimuthally modulated for an 
even number of electrons in a shell and round for an odd 
number per shell. The onset of this modulation is en- 
hanced within HF, which leads to overestimate the value 
of the critical density for the crossover to the Wigner 
molecule. However the HF densities mirror the classi- 
cal filling scheme with the electrons arranged in spatial 
shells. 
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TABLE I. Lowest energies of exact diagonalization and 
HF for N = 2 and various Sz and coupling strengths A. Rel- 

in %. 
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FIG. 1. N = 5. Absolute energy differences to QMC 
ground state (5* = 1/2), AE = Es — Eq^^ , for various spins 
vs. coupling constant A 




-4 -3 -2 -1 1 2 3 4 -4 -3 -2 -1 1 2 3 4 




-4 -3 -2 -1 1 2 3 4 -4 -3 -2 -1 1 2 3 4 
FIG. 2. Shadowed contour-plot of the HF density n^^ for 



\ = 6, Sz = N/2 and different electron numbers. Contours 
lie at integral multiples of 0.1 times the maximal density, (a) 
N = 2,{h)N = 3, (c) N = 4,{d)N = 5. 
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FIG. 3. 7-electron Wigner molecule. HF density n^^ for 
(a) \ = 2 S, = 1/2, (b) \ = 10 S,= 7/2. 
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FIG. 4. N = 6 Sz = 0. Rearrangement from a 6-fold 



[(a) A = 2] to a (deformed) 5-fold geometry [(c) A = 4, 
E = 41.509, (d) A = 6]. In (b) 6-fold isomer at A = 4 with 
energy E* = 41.838. 
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TABLE II. Energies of QMC vs. HE for iV = 8. Brack- 
eted numbers denote statistical errors of QMC. Relative error 

(i;UHF_^QMC)/^QMC 

A gQ^^ g"^ rel.err. 

2 4 48.3(2) 48.534 0.5 

2 3 47.4(3) 48.336 2.0 

2 2 46.9(3) 48.243 2.9 

2 1 46.5(2) 48.132 3.5 

4 4 69.2(1) 69?735 OS" 

4 3 68.5(2) 69.783 1.9 

4 2 68.3(2) 69.826 2.2 

6 4 86.92(6) 87.957 1.2 

6 3 86.82(5) 87.999 1.4 

6 2 86.74(4) 88.039 1.5 

8 4 103.26(5) 104.492 1.2 

8 3 103.19(4) 104.520 1.3 

8 2 103.08(4) 104.547 1.4 
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